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We consider 2-state quantum walks (QWs) on the line, which arc defined by two matrices. 
One of the matrices operates the walk at only half-time. In the usual QWs, localization 
does not occur at all. However, our walk can be localized around the origin. In this 
paper, we present two limit theorems, that is, one is a stationary distribution and the 
other is a convergence theorem in distribution. 

Keywords: limit distribution; localization; 2-state quantum walk. 



1. Introduction 

The 2-state quantum walk (QW) on the line Z = ( ] ±1, ±2, ■ ■ ■ } has been in- 
tensively studied, and the limit theorems are obtainedP^EllE] p or example, the 
limit distribution of the usual walks was calculatedP^ In the present paper, we 
consider a localization model of 2-state QWs. The motivation is the analysis of 
the time-inhomogeneous QWs. Our walks are determined by two matrices, one of 
which operates the walk at only half-time. The walk can be cons idered as one of the 
time-dependent models, for which there are some resultsP^ Particularly, Ref. E 
and ID discuss localization. We present the two limit theorems that show the local- 
ization of the probability distribution. One is calculation of the limit value for the 
probability P(X t — x) which walker is at position x £ Z starting from the origin, 
where t G {0, 1,2,.. .} is time, and the other is the convergence in distribution. In 
the usual walks, localization can not occurs. However, if we change the matrix at 
only half-time, then we find that the localization occurs from the results in this 
paper. The loca lization of QWs. w hich can be applied to quantum searchp] is of- 
ten investigated! I I I I I ' If limsup^^ P(X t — x) > for a position x, we 
call that the localization occurs. Therefore, our result insists that the localization 
occurs for any initial state. Moreover, we obtain the convergence in distribution of 
Xt/t as t — > oo. This limit distribution is described by both S- function and a den- 
sity function. For 3-state Grover walk, similar limit theorems were shownff^l The 
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limit distribution o f a 4-s tate walk corresponding to the 2-state walk with memory 
was also computed . SSH21 The present paper is organized as follows. In Section 2, 
we define our walk. We present the limit theorems as our main result in Section 3. 
Section 4 is devoted to the proofs of the theorems. By using the Fourier analysis, 
we obtain the limit distribution. Summary is given in the final section. 



2. Definition of a localization model of 2-state QWs 

In this section, we define a localization model of 2-state QWs on the line. Let 
| at) (x G Z) be an infinite components vector which denotes the position of the 
walker. Here, x-th component of |x) is 1 and the other is 0. Let \ipt(x)) G C 2 be the 
amplitude of the walker at position x at time t. The walk at time t is expressed by 

|* t )=5>)®|V> t (*)>. (1) 
The time evolution of our walk is depicted with the following two unitary matrices: 

(2) 

(3) 



u = 


cos 9 sin 9 
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s 














) 




sin 9 — cos 9 






•S 


— c 




H = 


cos 9\ sin 9\ 








Cl 


Sl 




sin 9\ — cos 9 


l _ 











G [0,2%), e ^ 0, f,7r, %E) and ci = cos6>i,si 



where c = cos 9, s = sm( 
sin#i (9\ € [0, 2tt)). Moreover, we introduce four matrices: 

P 



c s 



, Q = 


'0 " 

s — c 


, Pi = 


Cl Si 




,Qi = 


' 

Si -Ci 



Then, the evolution is determined by 

P\iH(x + l))+Q\iH(x-^)) (*^r) 
Pi \ijjt{x + 1)) + Qi \i/ft(x - 1)) (t = t) 



\i>t+i(x)) 



(4) 



(5) 



where r G {1, 2, . . .}. Note that P + Q = U and P x + Qi = H. The probability that 
the quantum walker X t is at position x at time t, P(X t = x), is defined by 



P(x t = x ) - (Mx)\Mx)) 



(6) 



In our main results, we focus on the probability distribution at time 2r + 1, 2r + 2. 
So, time r is called half-time in our walk. 

The Fourier transform |\fr t (fe)) (fc G [— 7T,7t)) of IV^Oe)) is given by 

\Mk)) = J2 e ~ ikX \^- ( ? ) 



By the inverse Fourier transform, we have 



liM*)> 



!**(*)> 



2tt' 



(8) 
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From §5§ and (J7J, the time evolution of \^ t (k)) becomes 

r u(k) i* t (fc)} (t + t) 

U(A)|* t (* ) (t = r) 



e Jfc 
e _i * 



where J7(fc) = R(k)U, H(k) = R{k)H and R(k) = 
that 

|*2r + l(fc)) =U(k) T H(k)U(k) T |* (fc)) , 
|*2r +2 (fc)) =U(k) T+1 H(k)U(k) T |# (fc)) 

In the present paper, we take the initial state as 

(T {a,f3] (a; = 0) 



lito(*)> = 



[0,0] (i^0) 



(9) 

Particularly, we see 

(10) 
(11) 

(12) 



where |a| 2 + \/3\ 2 = 1 and T is the transposed operator. We should note that 
\i&o(k)) = \ipo(0)). Figures [T]and [2] depict probability distributions of the walk under 
the condition 6 = ir/4, B\ = 0. In Figure [T] (a) , the probabilities at the position 
x = 0, ±2 are higher than other positions. Figure [2] shows the time evolution of 
probability distribution as r = 24. The walk evolves with the matrix H at only 
time t. So, the walk is a usual QW till time r. 

In comparison with both Figures [1] and [21 we show the probability distributions 
of a usual QW with 6 = 6\ = n/4 in Figures [3] and S) 
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i) |^ (0))= T [1/V2,VV2] 



(b) hMo)) = T [i,o] 



Fig. 1. The probability distributions at time t = 500 as r = 249 and 9 = 7r/4, $i = 0. 
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(b) |Vo(0))= T [l,0] 



Fig. 2. The evolution of the probability distributions for time t by density plot as r = 24 and 
6 = 7r/4,0i = 0. 
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Fig. 3. The probability distributions at time t = 500 with = 9i = ir/4 
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(a) \M0))= T [l/V2^/V2] (b) hM0))= T [l,0] 

Fig. 4. The evolution of the probability distributions for time t by density plot with 8 = 0\ = it/4. 



I 



I 



June 23, 2011 0:9 PRINTED 



Localization Model of 2- State QWs 5 



3. Limit theorems for the walk 

In this section, we show our main results. For our localization model of 2-state 
QWs, we obtain the two following limit theorems. 

Theorem 1. 

(i) For odd time 2t + 1. 



f {ciS - Sl S {1 ~ ]sir Ki(-l;P,a) (x = -l), 
(cifl - 81c e ) B 2(1 - |fl|)2 jCi(l;a,/3) (x = l), 



lim P(X 2t +i = x) = < 

T— >00 



2(c lS ~ SlC y 



]± ^My ]x] K 2 {x-,p,a) (x = -3,-5,-7,...), 



c*(i-M) 

2 %-^ S (^) 2,3:1 ^ «, P) (* = 3, 5, 7, . . .), 
(z = 0,±2,±4, 



(13) 



where 



#i(a;;a,/3) =c 2 \a\ 2 + 2s 2 (l - \s\)\(3\ 2 + sign(x)cs(l ~ \s\)(af3 + af3), (14) 
K 2 (x;a,(3) =c 2 s\a\ 2 + s(l - |s|) 2 |/3| 2 - si ff n(a;)c|s|(l - \s\){a0 + a/3), (15) 

and sign(x) means a sign of x. 
(ii) For even time 2t + 2, 

(£i£z£i^H)! ^3(_2; 0, «) (x = -2), 
(ci*-»iff (i-H) 2 ^. aj /3) (x = 2)) 

2|x| 



lim P(X 2r+ 2 = .x) = < 

T — ^OO 



2(e 1 s-s 1 c)^|s| 
7i 



2(c 1 s-s 1 c) 2 \s 



X 4 (x;Aa)(i = -4,-6,-8,...), 
(id^ 2|a V 4 ( x;aj/3 ) (x = 4,6,8,...), 
(a; = ±l,±3,...), 



(16) 
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K 3 (x;a,(3) =c 2 {l - S 2 | S |(2 - | s |)} \a\ 2 + 2s 2 (l - \s\) 3 \(3\ 2 

+ sign(x)cs(l + s 2 ){l ~ \s\) 2 (ap + a/3) , (17) 

K 4 (x;a,/3) =s 2 (|a| 2 - \/3\ 2 ) - sign{x)cs{a]3 + a/3) + \s\. (18) 



As a simple example of the case when the localization occurs in our walk, if 9 
tt/4,01 = and |V>o(0)) = T [l/V2 ,i/V2], then we have 



13 — 9\/2 

lim P(X 2 r+i = ±1) = = 0.136039 • • • 

O O fey 

lim P(X 2r+2 = 0) =— — — = 0.0857864 • • • 

T— fOO 2 

lim P(X 2T+2 = ±2) = 139 ~ 98 ^ - 0.101768 • 



(19) 
(20) 
(21) 



Figure [5] corresponds to the behavior of each probability in (fl~9)) , (|20|) and ([2Tj) . 

Next, we present the theorem of the convergence in distribution for X t /t, where 
t € {2r + 1, 2t + 2}. Some similar result s corresponding to Theorem 2 were shown 

[To 1 o 

for a 3-state walk or a 4-state walk! 1 1 Moreover, localization of multi-state walks 
was computed P^ The limit distribution of the usual 2-state walk does not have the 
delta measureP^ However, we find that the limit distribution of the 2-state walk 
defined in this paper has a delta- measure from Theorem 2. 
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Fig. 5. The behaviors of the probability P(X t = 0, ±1, ±2) as 9 = tt/4, 0i = 0. 
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Theorem 2. For our localization model of 2- state QWs, we have 

lim P(X t /t < x) = [ f(x)dx, (22) 

where 

f{x) =A5 (x) + 



7r(l — x 2 )\/c 2 — X 2 

I(-\c\,\c\)(x), 



l-<\a\ 2 -\8\ 2 



a 2 .x 4 + a\x 2 + a 



(23) 



c 2 (l-x 2 ) 

and So(x) denotes the delta-measure at the origin and Ia(x) = 1 if x € A, Ia{x) = 
if x A. The values A, ao,a\,a2 are independent on initial state |^o(0)) as follows: 

(cis - sic) 2 



A 



1 + \s\ 



(24) 



a =c . 



ai =2s 1 c(cis — sic) — c\, 



(25) 
(26) 

a 2 =(cis - sic) 2 . (27) 

Particularly, if 9\ = 0, then we obtain the density function of the usual walk, that 
is, 



/(*) 



7r(l — x 2 )\/c 2 — X 2 



1 



|2 (ap + a/3)s 



I(-\ c \,\c\)(x). (28) 



Figure E] shows the density function f(x) with 8 — tt/A,9i — 0. We have A = 
in Figure H 

We should note that the following relation between Theorems 1 and 2: 

V lim P(X t = x) = A. (29) 

Therefore, lim t _). 00 P(X t = x) is not a probability measure. 



-0.5 



(a) |^ (0)) = T [1/V2,z/V2] 



-1 -0.5 



(b) kMo)) = r [i,o] 



Fig. 6. The limit density function as 8 = 7r/4, Si = 0. 
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4. Proofs of theorems 



In this section, we will prove Theorems 1 and 2 in Section 3. Our approach is based 
on the Fourier analysis!^ 

4.1. Proof of Theorem 1 

At first, the eigenvalues Xj(k) (j — 1, 2) of U(k) can be computed as 



Ai(fc) = V 1 — c 2 sin 2 k + icsinfc, \2(k) = -yl — c 2 sin 2 k + ics'mk. (30) 
The normalized eigenvector \vj(k)} corresponding to \j(k) is 



\vj(k)) 



\/l — c 2 sin 2 fe ± c cos fc 
^ 2s 2 v / l-c 2 sin 2 /fc 



± a/ 1 — c 2 sin 2 — c cos /c 



(31) 



Therefore, the Fourier transform \^>o(k)} is expressed by \vj(k)) as follows: 



\i>o(k)}^J2^( k ^o(k))\v J (k)). 

3=1 



(32) 



In the proof, we focus on even time 2r+2. From (fTTj) and (|3"2"1) . the Fourier transform 
at time 2r + 2 is given by 

|*2r +2 (fc)> =u{k) T + 1 H{k)u{ky |* (*)) 

2 

=J7(fc) T+1 ff(fc) ^ AJ t (fc) (^(fc)|* (fc)> |^(fc)) 

3'i=l 

2 

= E W («ii(*)l*o(*)> (j/(fc) T+1 i?(fc) \v h (k))) ■ (33) 

3'i=l 

Moreover, rewriting U(k) T+1 H(k) \vj 1 (k)) as 

2 

U(ky +1 H(k) \v n (k)) =U(k) T+1 J2 (v 32 (k)\H(k)\v n (k)) \v j2 (k)) 

32 = 1 

2 

= E A J, +1 ( fc ) («A(*)l^(*)l«ii(*)> M fc )> . (34) 

32=1 



we obtain 



|*2r +2 (fc)> = E E ^l(*) T ^(*) r+1 («ii(*)l*o(*0) 
31=1 32=1 

x (^(fc)!^)!^^))!^^)). 



(35) 
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We should note \x(k)\2(k) = — 1. Calculating the inverse Fourier transform, we 
have 



2 2 

\4>2r+2(x)) =E E T X nW\ j2 (kY +1 (V n (k)\* (k)) 
.71=1.72 = 1 J ~ lt 



Akx dk 



x (v h (k)\H(k)\v h {k)) \v h (k)) e lfc * — . (36) 



By using the Riemann-Lebesgue lemmaj^ we see 



|^2r+2(x)) ~ (-1) T / Al(fc) (v 2 (fc)|*o(*0) ( Vl (fc)|ff(fc)| U2 (fc)) e 



//,•.(• 



(-l) r / A 2 (fc) («i(fe)|* (fc)> M*)l# (*)M*)> M*0> e 



2^ 

2tt' 
(37) 



where g(r) ~ /i(r) denotes lim T _ i . 00 g(r)//i(r) = 1. From ([37|). we get 



\^2t+2{x)) ~ < 



(-i)T( eiS -s lC )H(i-H) 
^ 



(-l) T (cis-sic) 



(-l) T (ciS-Oie) 



^ 



-csl 2 a-(l-\s\)l3 
I 2 {\s\(l-\s\)a + csl3} 

I 2 {csa-\s\(l-\s\)(3} 
(1 - |s|)q - csl 2 p 

±csa - |s|(l =F \s\)P~ 
\s\(l± \s\)aTcsf3 _ 

"0" 





i = 0), 
i = -2), 
1 = 2), 
x = ±4, ±6 



as = ±1,±3 



(38) 
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where I x = { ^j^p^ | ■ Similarly we can compute \ip2T+i(%)) a s follows: 



|^2r+i(a0) 



(-iy( cis - si c)(i-\ s \) 
3 



(-iy( Cl s- Sl c)(i-\ s \) 

7S. 



c(sa — cj3) 
-|s|(l- \s\)a-csfi 

csa- \s\(l- \s\)0 
—c(ca + s/3) 



-c 2 {s(l- \s\)a + c\s\/3} 
(l-\s\){c\s\(l-\s\)a + c 2 s0} 

(l-| s |){-c2 sa + c| S |(l-| S |)^} 
-c 2 {c| S |a-.s(l-| S |)/3} 



(x = -1), 
(x = 1), 

(x = -3, -5, . . .), 

(x = 3,5, . . .), 

(.x = 0,±2,±4,...), 
(39) 



where J x = ^2T^r7xzq7rf ■ From (|3"5)) and (|3"5)) . the proof is completed. 



□ 



4.2. Proof of Theorem 2 

We calculate the characteristic function E(e tzXt ^ t ) as t — > oo, where E(X) denotes 
the expected value of X. At first, (|35p can be written as 



|*2r +2 (fc)) =u(ky +1 H(k)u(ky |* (fe)> 

=A 1 (fc) 2T+2 ^iW |^(fc)) + A 2 (fc) 2T+2 A 2 (fc) |«2(fc)) 

- A 3 (fc) 2r+2 A 3 (fc) | Wl (fc)) - A 4 (fc) 2r+2 A 4 (fc) h(fc)> , (40) 



where A 3 (fc) = A 4 (fc) = i and 

Ai(fc) =Ai(fc) ( Vl (fc)|^ (fc)) , (41) 

A 2 (fc) =A 2 (fc) ( U2 (fc)|# (fc)) («2(fc)|H(fc)|«2(A;)) , (42) 

A 3 (fc) =Ai(fc) (w 2 (fc)|#o(fc)) («i(fc)|#(fc)MA;)) , (43) 

A 4 (fc) =A 2 (fc) ( Vl (fc)|# (fc)) («2(fe)|^(fc)|«i(fc)> ■ (44) 
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Substituting 2r + 2 = t, we can calculate the r-th moment of X t as follows: 

E((x t y)=^ r P(x t = x) 



(* t (fc)i(D'-i*t(fc)r 

-7T 



2tt 



=(t) r /* J ^fy(fc) r |A,f - A 1 (fc)*A 3 (fc) t /i 1 (fc) r , Ai(fc)i4 3 (fc) 
- AT(fc) t A 3 (fc) t /i3(fc) r 3r(fc)A3(fc) - A 2 (fc) t A7(fc) t ^ 2 (fc) r A 2 (fc)3;(fc) 

^+0(t r - 1 ), (45) 

where ^(fc) = DXj(k)/ Xj(k), D = i(d/dk) and (t) r = t(t - 1) x • • • x (t - r + 1). 
By using the Riemann-Lebesgue lemma, we have 

^]T^(fc)|^(fc)| 2 - 

2 

=° rA + /^E^( fc )l^( fc )| 2 §- W 



where 

2 dfc (cis — sic) 



«/ — 7T • o 



2tt 1+ s 

J=3 



Therefore, we obtain 



lim £((X t /t) r ) =0 r A+ / x r ^ - 

t^°° J-oo 7r(l - X 2 )VC 2 - X 2 



1 - < hi 2 - I/3I 2 + rz_i_zi^ u 



(47) 



a/3 + S^)s ' 



a 2 x 4 + aix 2 + ao 



/(-| c |,| c |)(a;)rfa; 



c 2 (l-.x 2 ) 

/OO 
x r f(x)dx, (48) 
-oo 



where 



a =c 2 , (49) 
ai —2sic(cis — sic) — cl, (50) 
a 2 =(cis - sic) 2 . (51) 

We should remark 

/oo 
0(a;)5o(aOda; = s(O). (52) 
-OO 
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12 Takuya Machida 

By (|48|) . we can compute the characteristic function E{e lzXt ^ t ) as t — > oo. Thus the 
proof of Theorem 2 is completed. 

□ 

5. Summary 

In the final section, we conclude and discuss the probability distribution of our 
walks. In the usual 2-state walk defined by the matrix U , the localization does not 
occur at all. However, if another matrix H U) operates the walk at only half-time, 
then localization occurs. In Theorem 1, the behavior of the probability P(X t — x) 
was calculated as t — > oo. Moreover, we found that the limit distribution of X t /t 
had both a delta measure and a density function from Theorem 2. The interesting 
problem is calculation of the limit distribution for the walk which the matrix H 
operates more than twice. 
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